MitacsTSFCanadaLocations

library(tidyverse)

## ── Attaching packages ─────────────────────────────────────────────────────────────────────────── tidyverse 1.2.1 ──

## ✔ ggplot2 3.2.1     ✔ purrr   0.3.3
## ✔ tibble  2.1.3     ✔ dplyr   0.8.3
## ✔ tidyr   1.0.0     ✔ stringr 1.4.0
## ✔ readr   1.3.1     ✔ forcats 0.4.0

## ── Conflicts ────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()

library(maps)

## 
## Attaching package: 'maps'

## The following object is masked from 'package:purrr':
## 
##     map

library(maptools)

## Loading required package: sp

## Checking rgeos availability: TRUE

library(sp)
library(mapdata)
map('worldHires','Canada', xlim=c(-141,-53), ylim=c(40,85), col='gray90', fill=TRUE)

df = read.csv('https://raw.githubusercontent.com/bcdunbar/datasets/master/meteorites_subset.csv')

Sys.setenv('MAPBOX_TOKEN' = 'pk.eyJ1Ijoic2FsbHlpbm5pcyIsImEiOiJjazR3eWszdnEwN3B6M250aGt3YnM4aW9yIn0.KQelI9-fKQw9FfAkQRo6jA')

library(readxl)
CadTSF <- read_excel("Canadian Filtered TSFs.xlsx", 
    sheet = "Filtered and Changed")

## New names:
## * `` -> ...4

#View(CadTSF)

library(scales)  #for transparency

## 
## Attaching package: 'scales'

## The following object is masked from 'package:purrr':
## 
##     discard

## The following object is masked from 'package:readr':
## 
##     col_factor

map('worldHires','Canada', xlim=c(-141,-53), ylim=c(40,85), col='gray90', fill=TRUE)
points(CadTSF$Longitude, CadTSF$Latitude, pch=19, col="red", cex=0.5) 

library(mapproj)
map('worldHires','Canada', xlim=c(-141,-53), ylim=c(40,85), col="grey80", fill=TRUE, projection="gilbert", orientation= c(90,0,225))
lon <- c(-72, -66, -107, -154)  #fake longitude vector
lat <- c(81.7, 64.6, 68.3, 60)  #fake latitude vector
coord <- mapproject(CadTSF$Longitude, CadTSF$Latitude, proj="gilbert", orientation=c(90, 0, 225))  #convert points to projected lat/long
points(coord, pch=20, cex=.8, col="red")  #plot converted points

library(plotly)

## 
## Attaching package: 'plotly'

## The following object is masked from 'package:ggplot2':
## 
##     last_plot

## The following object is masked from 'package:stats':
## 
##     filter

## The following object is masked from 'package:graphics':
## 
##     layout

packageVersion('plotly')

## [1] '4.9.0'

library(plotly)

dat <- map_data("world","canada") %>% group_by(group)

p <- plot_mapbox(dat, x = ~long, y = ~lat) %>%
  layout(mapbox = list(zoom = 4,
                       center = list(lat = ~median(lat),
                                     lon = ~median(long))
  ))

p

## No scattermapbox mode specifed:
##   Setting the mode to markers
##   Read more about this attribute -> https://plot.ly/r/reference/#scatter-mode

library(plotly)
library(dplyr)

p <- plot_mapbox(mode = 'scattermapbox') %>%
  add_markers(
    data = CadTSF, x = ~Longitude, y = ~Latitude, text= ~Mine, pch = 3, color= ~Status, size = ~Current_Maximum_Height_m, hoverinfo = "text") %>%
 layout( plot_bgcolor = 'white', paper_bgcolor = 'white',
         mapbox = list(zoom = 3,
                       center = list(lat = 60,
                                     lon = -97)
  ))
p